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ABSTRACT 

Motivated by the large-scale asymmetry observed in the cosmic microwave background sky, we 
consider a specific class of anisotropic cosmological models - Bianchi type VIR - and compare them 
to the WMAP first-year data on large angular scales. Remarkably, we find evidence of a correlation 
which is ruled out as a chance alignment at the 3 a level. The best fit Bianchi model corresponds to 
x = 0.55, fl 0 = 0.5, a rotation axis in the direction (Z,6) = (222°,—62°), shear (jj) 0 = 2.4 x 10“ 10 
and a right-handed vorticity (j[) 0 = 6.1 x 10” 10 . Correcting for this component greatly reduces the 
significance of the large-scale power asymmetry, resolves several anomalies detected on large angular 
scales (ie. the low quadrupole amplitude and quadrupole/octopole planarity and alignment), and 
can account for a non-Gaussian “cold spot” on the sky. Despite the apparent inconsistency with the 
best-fit parameters required in inflationary models to account for the acoustic peaks, we consider the 
results sufficiently provocative to merit further consideration. 

Subject headings: cosmic microwave background — cosmology: observations — methods: numerical 


1. INTRODUCTION 

Despite initial reactions to the first year WMAP * 1 2 3 4 5 data 
- “the most revolutiona ry result is that there are no revo¬ 
lutionary results” llBahca,111200311 - continued assessment 
has revealed some interesting discrepancies with the in¬ 
ferred best-fit, cosmological constant dominated, infla¬ 
tionary model of the primordial fluctuation spectrum. 

A debate about the apparently low quadrupole am¬ 
plitude (ori ginall y observed by COBE-DMR) has arisen 
lEfsta,thionll2004D along with cla ims ab o ut other anoma¬ 
lies on large angular scales. Ide Oliveira-Costa et all 
(2004) demonstrated a curious planarity an d alignment 
of the quad rupole and octopol e. In addition. iVielva et al.l 
1200411 and I Cruz et aL (|20f)5f l have detected a localised 
source of non-Gaussianity, in the form of a very cold 
spot on the sky of angular scale ~ 10°. Of particu¬ 
lar interest to us, however, was the discovery of a re¬ 
markable asymmetry in large-scale power as measured 
in the two hemispher e s of a particular referen ce frame. 
lEriksen et all l2004al I Hansen et~a!1 12004 albll. It re¬ 
mains unclear why this was not detected b y the bipo¬ 
lar power spectrum of lHaiian. Souradeep. & Cornish 1 
1200511 . Nevertheless, this has provided the motivation 
to investigate anisotropic cosmo logical models. 

Following i Bunn. Ferreira, fc Silk lilOOfilf and 

iKogut et al . 1199711 . we foc us on a specific c lass of 
models - Bianchi type VII/* ( Bar row et alJH985h . Such 
models were previously compared to the COBE-DMR 
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data in IKogut et alJ : 199 71. where limits on the shear 
(jf) 0 < 10 -9 and vorticity (jf) 0 < 6 x 10~ 8 were estab¬ 
lished. We consider that the observed anisotropy is the 
sum of two contributions - an ‘isotropic’ term which is 
connected to variations in the density and gravitational 
potential, and a term from the anisotropic metric. A 
mechanism to generate the isotropic fluctuations is 
required, particularly on small angular scales where the 
Bianchi contribution is negligible, but this need not 
be inflationary. In our analysis we smooth the data to 
probe only the low-£ regime, and find that results are 
practically insensitive to the choice of spectrum for the 
isotropic component on large-angular scales. 

Remarkably, we find a statistically significant correla¬ 
tion between one of these models and the WMAP data. 
Such a result may help to resolve some of the more un¬ 
usual observed features of the microwave sky, albeit at 
the introduction of a new conundrum - the large-scale 
anisotropy is described at least in part by a low density 
Bianchi model, whereas the smaller scale fluctuations are 
consistent with a cosmological constant dominated, crit¬ 
ical density model in an inflationary scenario. 

2. DATA, SIMULATIONS AND TEMPLATES 

In this analysis, we utilize the first-year WMAP 
data 6 i n both raw and foregr ound template corrected 
forms iBennett et alJ l2f)03albH and three heavily pro¬ 
cessed maps derived to minimize the foreground con¬ 
tribution in a temp late in dependent wav (WI LC - 
|Bennett_e^lJl2003li LILC - lEriksen et aljl20fl4bt TOH 
UdhoTnariTTtaTl 1 2 D 03ll . For the three ILC-like maps we 
consider the f ull sky in the analysi s, whereas the Kp2 
or KpO masks ( IBennett et alJl2003bll are imposed other¬ 
wise. All are smoothed to a common resolution of 5.5° 
FWHM and downgraded to a HEALPix' resolution of 
A^ide = ^2 corresponding to harmonics up to £ = 64. 
Such processing does not compromise the analysis since 

6 available at http: //lambda.gsfc.nasa.gov 

7 http://www.eso.org/science/healpix/ 
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the Bianchi models exhibit little power outside of this 
Grange. 

To assess the statistical significance of the fits we 
compare our results to an ensemble of 10 000 LILC 
simul ations, produced using the pipeline described by 
lEriksen et all d2004bl) . 

The Bianchi mo del tem pl ates ar e constructed using the 
formalism of lBarrow et al . (jl985t) . These models are pa¬ 
rameterized by fio, x = — fi„) (where h is the 

scale on which basis vectors change orientation), and a 
handedness. These models have a prefered axis along 
which the expansion rate is different, and about which 
the basis vectors themselves rotate. The smaller the x, 
the more rotation in a given distance traversed along 
a geodesic, and therefore the tighter the observed spi¬ 
ral pattern of induced temperature anisotropies. The 
smaller the STq , the larger is the asymmetry along the axis 
and the more focused is the structure in only one direc¬ 
tion. In this analysis, we consider the range 0.1 < x < 10 
and 0.1 < S2o < 1. The resulting temperatu re anisotr opy 
patter n is described by Equation 4.11 of lBarrow et all 
( 1985 ). The fit amplitude determines the shear (jj) 0 
and the vorticity (jj) 0 following Equation 4.8. 

3. METHOD 

In order to compute the correlation betwe en the data d 
and a temp late t, we adopt the formalism of Gorsk i~et all 
(119961) and lKogut et all ; 1997 ) . Specifically, at a given 
frequency we estimate a template coupling constant a by 
minimizing y 2 = (d — at)^ M _1 (d — at), where M is the 
CMB signal plus noise covariance matrix. The solution 
to this problem (and corresponding uncertainty) is a = 
(t T M _1 d)/(t T M^ 1 1) and 6a = (t T M -1 1) -1 / 2 . The 
generali zatio n to the multi-frequency case can be found 
in IGorski et alJ (I1996D . In what follows, the CMB signal 
covariance is specified by the WMAP best-fit theoretical 
power spectrum. 

In our analysis, we must determine a over all possible 
models and, for each model, all relative rotations between 
the corresponding template and the data. Such an anal¬ 
ysis is generally computationally demanding, but can be 
accelerated by a full-sky harmonic-space approach. In 
this case, if the noise covariance is a diagonal approxi¬ 
mation based on the mean noise level for the map under 
consideration 8 then the denominator of a is invariant 
under template rotation and simply computed for any 
given model. The numerator can then be evaluated from 
the convolution of the template and the “whitened” data 
map M~ 1 d. This oper ation is perform ed usi ng th e total 
convolution method of lWandelt and Gorskil 6 2001 5 to ef¬ 
ficiently evaluate all such convolutions over all possible 
rotations. 

The technique is only straightforward when consider¬ 
ing full sky coverage. Thus, despite concerns about the 
effectiveness of foreground removal in the Galactic plane, 
we analyse the foreground-cleaned LILC map. Since 
chance alignments between the Bianchi template and the 
CMB anisotropy arising even in the absence of rotation 
and shea r ca n gene rate non-zero correlations, we follow 
(Kogut et al . lfl997t) and define the parameter T = a/5a 

8 This approximation has a negligible effect on results at the 
WMAP signal-to-noise ratio. 


to help reject this possibility at a given confidence and 
to select the best model. A reference ensemble of 10 000 
LILC simulations were analysed using the above formal¬ 
ism to establish the distribution of a values and hence 
to estimate the significance of our best fit model. 

Once the model and orientation axes are specified, we 
apply a pixel-based analysis to various combinations of 
the WMAP sky maps and impose the conservative KpO 
Galactic cut to reduce foreground contamination, in or¬ 
der to confirm these results. In the case of incomplete 
sky analysis, the amount of template structure which is 
masked can affect the uncertainty of the fit. Neverthe¬ 
less, for the best fit orientation where the structure is 
largely outside the excluded Galactic Plane region, we 
find that the fit and its uncertainty are dominated by 
the CMB covariance rather than the particulars of the 
noise properties or sky cut. Therefore, although we cal¬ 
ibrate using simulations strictly only applicable to the 
LILC result, we expect that the predicted levels of chance 
alignment detections should be quite robust. 

4. RESULTS 

The results from the analysis are summarized in Ta¬ 
ble U for our best fit model at x = 0.55 and fl 0 = 0.5, 
right-handed. The excess power in the southern hemi¬ 
sphere requires a model with low Hq to focus the spi¬ 
ral. The cold spot and surrounding structure effec¬ 
tively select the particular x. See Figure^ The best- 
fit Euler angles 9 for the WILC map were found to be 
($ 2 , 0, 4>i) = (42°, 28°,—50°), and in what follows this 
is adopted as the best-fit axis. However, both the LILC 
and TOH sky maps as analysed on the full-sky yielded 
statistically consistent best-fit axes. This rotation places 
the center of the spiral structure, originally at the —z 
axis, at (/,&) = (222°,—62). For each map, we tabulate 
the sheer and vorticity, and corresponding estimate of 
significance as computed from the fraction of LILC sim¬ 
ulations with a smaller a value than the observed map. 
The estimated shear amplitudes are consistent with a 
value ~ 2.4 x 10~ 10 , significant at the 99.8% confidence 
level relative to the Monte Carlo ensemble. 

In order to assess the effect of residual foreground con¬ 
tamination (ie. due to signals not well-traced by the 
foreground templates) several difference maps were con¬ 
structed from the WMAP foreground corrected maps. 
The results indicate that foreground residuals are un¬ 
likely to result in spurious detections. 

In Figure |T| we show (from top to bottom) the WILC 
map, the best-fit Bianchi model, and the difference be¬ 
tween the two. It should be apparent that the “Bianchi 
corrected” map exhibits greater isotropy than the WILC 
data. 

In Figure |2 we compare the power spectra of the orig¬ 
inal and the Bianchi-corrected V+W linear combina- 
tion map 10 , as co mputed by the MASTER algorithm 
llHivon et alil2002 ) for the Kp2 sky coverage. It should 
be noted that, whilst the quadrupole amplitude is signif¬ 
icantly increased, on average the Bianchi-corrected map 
has slightly less power on large scales, although both 
corrected and uncorrected data are in good agreement 

9 We adopt the zyz-convention of lWandelt and Odrski ( 200111 for 
the three Euler angles. 

10 The WMAP team <Hinsha,w et a,l.I200H1 use data solely from 
the V- and W-bands over the spectral range t < 100. 
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Fig. 1. — Top panel: WMAP Internal Linear Combination map. 
Middle panel: Best-fit Bianchi VIIjj template (enhanced by a factor 
of four to show the structure). Bottom panel: Difference, i.e. the 
“Bianchi-corrected” ILC map 


TABLE 1 

Fitted template amplitudes 


Map 

(i) 0 

(xio -10 ) 

(#)o 

(Xl0 _r °) 

^(“sim < a obs) a 

WILC 

2.39+0.45 

6.14 

99.7% 

LILC 

2.37+0.45 

6.07 

99.6% 

TOH 

2.25+0.45 

5.76 

98.6% 

K b 

2.39+0.46 

6.14 

99.7% 

Ka b 

2.35+0.46 

6.04 

99.5% 

Q b 

2.35+0.46 

6.04 

99.5% 

V b 

2.42+0.46 

6.22 

99.8% 

W b 

2.48+0.46 

6.38 

99.9% 

Q+V+W bc 

2.42+0.46 

6.22 

99.8% 

V+W bc 

2.43+0.46 

6.25 

99.8% 

Q-V bc 

-0.07 + 0.01 

0.17 

- 

v_w bc 

-0.06 + 0.01 

0.14 

- 

Q_W bc 

-0.12 + 0.01 

0.31 

- 


Note. — Amplitudes of the best fit Bianchi model. See 
text. a Fraction of LILC simulations with lower amplitudes. 
b Including simultaneously derived foreground corrections. c 
Using WMAP frequency maps pre-corrected for Galactic fore¬ 
grounds. 


for t > 15. More importantly, the general shape of the 
Bianchi-corrected spectrum is flatter than for the WMAP 
best-fit power spectrum, and i n remark able agreem ent 
with the theoretical fit made hv lHansen et a.U ( 2004al ) to 
the northern hemisphere data (defined in the reference 
frame which maximises the power asymmetry) which 
favours a lower value for r. 

5. IMPLICATIONS 



0 10 20 30 

Multipole / 


Fig. 2. — Comparison of power spectra. The gray and black solid 
lines show the power spectrum estimated from the co-added V+W 
map before and after correcting for the Bianchi template, respec¬ 
tively. The dotted gray and black lines shows the th eoretical best- 
fit powe r-spectra from the WMAP -team analysis and IHansen et alj 
(: 2004a ) respectively. The latter is a fit to northern hemisphere data 
alone. The light gray dashed line corresponds to the full sky power 
spectrum of the best-fit Bianchi template. 


Quadrupole amplitude: The WMAP team suggest 
that the quadrupole amplitude is significantly low, al¬ 
though other analyses have found it to be quite accept¬ 
able (Slosar & Seljak 2004; O’Dwyer et al. 2004). As 
seen in Figure O the Bianchi-corrected V+W map has a 
quadrupole amplitude of 504 pK 2 , compared to the un¬ 
corrected amplitude of 137 pK 2 and the WMAP theoret¬ 
ical best-fit spectrum value of 869 pK 2 . In this context, 
the quadrupole amplitude should no longer be consid¬ 
ered anomalous. Whether the amplitude enhancement 
itself requires an unusual cancelation between the intrin¬ 
sic and Bianchi-induced quadrupoles, which could also 
be considered a ‘fine-tuning’of the model, is deferred to 
a later analysis. 

Low-f anomalies: An alignment between the 

qua drupole and the octo p ole h as bee n claimed 
by jde Oliveira-Costa et al . li20ril and iCopi et all 
<200+ Quantitative ca lc ulatio ns sim ilar to those 
oFJdj^Ohveira-Costa et Id] <t200l and lEriksen et alJ 
ll2004bh show that a stronger planarity is expected by 
chance with a probability of 52% for both the £ = 2 and 
3 modes after subtracting the Bianchi template. The 
angle between the preferred directions of the 1 = 2 and 
3 modes is 70° after subtracting the Bianchi template, 
compared to 12° before. Additionally, the £ = 5 and 6 
modes (Eriksen et al. 2004b) become less anomalous, 
with significances dropping to 1.5tr in bot h cases. 
Large-scale power asymmetry: lEriksen et alJ 

: 200 1 al : reported that the large-scale power (£ < 40) in 
the WMAP data is anisotropically distributed over two 
opposing hemispheres, with a significance of 3er com¬ 
pared with simulations. Repeating the analysis and 
adopting the Kp2 sky coverage, we compare the cor¬ 
rected V+W WMAP map with 2048 simulations. We 
find that 13.6% of the simulations have a larger maxi¬ 
mum power asymmetry ratio than the Bianchi-corrected 
map, whereas only 0.7% have a larger ratio than the 
uncorrected data. It is apparent that the maximum 
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Wavelet scale (deg) 

Fig. 3. — The kurtosis of the wavelet coefficients as a func¬ 
tion of wavelet scale. The solid (dotted) lines show the results 
for the northern (southern) Galactic hemisphere (excluding the re¬ 
gion where \b\ < 20°) computed from the WMAP ILC map. The 
squares (crosses) show the kurtosis before (after) correction for the 
Bianchi template. Grey bands indicate the two-sided 68% and 95% 
confidence intervals obtained from 1000 Monte-Carlo simulations. 


power ratio between any two hemispheres is significantly 
suppressed after subtracting the Bianchi template - no 
asymmetry axis is found at any statistically significant 

level. _ _ 

Wavelet kurtosis: lYielva ~et"aLl ll200l used a wavelet 
technique to detect an unusually cold spot (~ 3<r sig¬ 
nificance relative to Gaussian simulations) at Galactic 
coord inates (l,b) = (207°, — 59°). We repeat the analy¬ 
sis o f lVielva et, al.l ( 2004 ). and compute the kurtosis of 
the wavelet coefficents as a function of scale from both 
the WILC and the corresponding Bianchi-corrected map. 
For computational convenience, a Galactic cut is imposed 
to exclude the region where \b\ < 20°. The results from 
this exercise are reported in Figure 0 After subtracting 
the Bianchi template, the southern hemisphere results 
are consistent with the 95% confidence intervals obtained 
from simulations - no non-Gaussian features are appar¬ 
ent. 

6. CONCLUSIONS 

We have considered the Bianchi type VII/j class of 
anisotropic cosmological models and fitted the predicted 
CMB temperature anisotropy patterns to the first year 
WMAP data. The results are essentially independent 


of frequency or the method of foreground subtraction. 
A particular model with a rotation axis in the direction 
(Z, 6) = (222°, -62°), shear (f) Q = 2.4 x 1(T 10 and a 
right-handed vorticity (jj) Q = 6.1 x 10 -10 yields the 
best fit to the data. An ongoing analysis will improve 
the accuracy and completeness of the search of the model 
space and examine other possibly significant models. 

Should this result be considered more than a curios¬ 
ity? A skeptic would no doubt stress the inconsis¬ 
tency between those cosmological parameters describ¬ 
ing the Bianchi anisotropy pattern (an open model with 
n 0 = 0.5) and the WMAP best-fit cosmological power- 
spectrum which must account for the acoustic peaks. 
Yet paradoxically, it is precisely the low value for Do 
which is required to allow the focussing of the Bianchi 
anisotropy pattern into one hemisphere. The preferred 
Bianchi model reconciles the observed asymmetric distri¬ 
bution of power on large angular scales, disrupts the ob¬ 
served planarity and alignment of the quadrupole and oc- 
topole moments, and provides an explanation for a highly 
non-Gaussian signature on the sky - unexpected results 
given that the model was not selected based on these 
criteria, but solely on a statistical fit of the predicted 
anisotropy pattern to the WMAP data. Furthermore, 
the model is a l so cons istent with the tentative result from 
lHansen et all ll2004aD that the estim ated optical depth 
of t = 0.17 on the (nearly) full sky llKognt et al]l20tbl 
could in large part originate in structure associated with 
the southern hemisphere (in the reference frame which 
maximizes the power asymmetry). 

While the consistency of our result on large angular 
scales with models for the small-scale anisotropy remains 
ambiguous, we note that the models considered include 
no dark energy contribution. A more self-consistent ap¬ 
proach is clearly warranted, but we propose that the re¬ 
sult is sufficiently provocative to encourage a renewed 
exploration of anisotropic models in a broader cosmolog¬ 
ical parameter space. 
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Data Analysis (LAMBDA). 


REFERENCES 


Bahcall, J. 2003, http://www.sns.ias.edu/~jnb/Papers/Popular/ 
MAPremarks.html 

Barrow, J. D., Juszkiewicz, R., & Sonoda, D. H. 1985, MNRAS, 
213, 917 

Bennett, C. L. et al. 2003a, ApJS, 148, 1 

Bennett, C. L. et al. 2003b, ApJS, 148, 97 

Bunn, E. F., Ferreira, P. G., Silk, J., 1996, PRL, 77, 2883 

Copi, C. J., Huterer, D., & Starkman, G. D. 2004, Phys. Rev. D, 
70, 043515 

Cruz, M., Martinez-Gonzalez, E., Vielva, P. &; Cayon, L. 2005, 
MNRAS, 356, 29 

de Oliveira-Costa, A., Tegmark, M., Zaldarriaga, M., Sz Hamilton, 
A. 2004, Phys. Rev. D, 69, 063516 

Efstathiou, G. 2004, MNRAS, 348, 885 

Eriksen, H. K., Hansen, F. K., Banday, A. J., Gorski, K. M., &z 
Lilje, P. B. 2004a, ApJ, 605, 14 


Eriksen, H. K., Banday, A. J., Gorski, K. M., $>z Lilje, P. B. 2004b, 
ApJ, 612, 633 

Gorski, K. M., Banday, A. J., Bennett, C. L., Hinshaw, G., Kogut, 
A., Smoot, G. F., & Wright, E. L. 1996, ApJ, 464, Lll 

Gorski, K. M., Hivon, E., Banday A. J., Wandelt, B. D., Hansen, 
F. K., Reinecke, R., & Bartelmann, M., 2005, ApJ, 622, 759 

Hansen, F. K., Balbi, A., Banday, A. J., &z Gorski, K. M. 2004a, 
MNRAS, 354, 905 

Hansen, F. K., Banday, A. J., & Gorski, K. M. 2004b, MNRAS, 
354, 641 

Hinshaw, G. F. et al. 2003, ApJS, 148, 135 

Hivon, E., Gorski, K. M., Netterfield, C. B., Grill, B. P., Prunet, 
S., & Hansen, F. 2002, ApJ, 567, 2 

Kogut, A., Hinshaw, G., Banday, A. J. 1997, Phys. Rev. D, 55, 4, 
1901 

Kogut, A., et al. 2003, ApJS, 148, 161 




















5 


O’Dwyer, I. J., et al. 2004, ApJ, 617, L99 
Slosar, A., Sz Seljak, U. 2004, Phys. Rev. D70, 083002 
Hajian, A., Souradeep, T. ,Sz Cornish, N., 2005, ApJ, 618, L63 
Tegmark, M., de Oliveira-Costa, A., Sz Hamilton, A. 2003, 
Phys. Rev. D, 68, idl23523 


Vielva, P., Martinez-Gonzalez, E., Barreiro, R. B., Sanz, J. L., Sz 
Cayon, L. 2004, ApJ, 609, 22 

Wandelt, B. D. & Gorski, K. M. 2001, Phys. Rev. D, 63, 123002 



